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Abstract 

We  investigate  analytically  and  numerically  the  effects  of  tailoring  the  density  profile  in  a 
laser  target  in  order  to  decrease  imprinting  of  mass  perturbations  due  to  the  long-wavelength  modes. 
Inverting  the  acceleration  of  the  ablation  front  during  the  shock  transit  time  could  reduce  the  early- 
time  mass  perturbation  amplitudes  developed  in  the  target  after  the  shock  transit.  This  principle  was 
first  suggested  for  mitigating  the  RT  instability  of  imploding  Z-pinches  [Velikovich  et  al.,  Phys. 
Rev.  Lett.  77,  853  (1996);  Phys.  Plasmas  5,  3377  (1998)].  As  the  shock  wave  slows  down 
propagating  into  higher  density  layers,  the  effective  gravity  near  the  ablation  front  has  the  same 
direction  as  the  density  gradient.  This  makes  the  mass  perturbations  near  it  oscillate  at  a  higher 
frequency  and  at  a  lower  amplitude  than  they  normally  would  due  to  the  “rocket  effect”  caused  by 
mass  ablation  [Sanz,  Phys.  Rev.  Lett.  73,  2700  (1994);  Piriz  et  al.,  Phys.  Plasmas  4,  1117 
(1997)].  So,  tailoring  density  profiles  instead  of  using  flat  densities  is  demonstrated  to  reduce  the 
“seed”  mass  perturbation  amplitude  at  the  onset  of  the  exponential  RT  growth. 

PACS  numbers:  52.40.Nk,  52.35.Py 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

1999 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-1999  to  00-00-1999 

4.  TITLE  AND  SUBTITLE 

Reduction  of  early-time  perturbation  growth  in  ablatively  driven  laser 
targets  using  tailored  density  profiles 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,LCP&FD, 4555  Overlook  Avenue 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

his  article  appeared  in  Physics  of  Plasmas  and  may  be  found  at  Phys.  Plasmas  6,  3283  (1999) 


14.  ABSTRACT 

We  investigate  analytically  and  numerically  the  effects  of  tailoring  the  density  profile  in  a  laser  target  in 
order  to  decrease  imprinting  of  mass  perturbations  due  to  the  long-wavelength  modes.  Inverting  the 
acceleration  of  the  ablation  front  during  the  shock  transit  time  could  reduce  the  earlytime  mass 
perturbation  amplitudes  developed  in  the  target  after  the  shock  transit.  This  principle  was  first  suggested 
for  mitigating  the  RT  instability  of  imploding  Z-pinches  [Velikovich  et  al.,  Phys.  Rev.  Lett.  77,  853  (1996); 
Phys.  Plasmas  5,  3377  (1998)].  As  the  shock  wave  slows  down  propagating  into  higher  density  layers,  the 
effective  gravity  near  the  ablation  front  has  the  same  direction  as  the  density  gradient.  This  makes  the 
mass  perturbations  near  it  oscillate  at  a  higher  frequency  and  at  a  lower  amplitude  than  they  normally 
would  due  to  the  "rocket  effect"  caused  by  mass  ablation  [Sanz,  Phys.  Rev.  Lett.  73,  2700  (1994);  Piriz  et 
al.,  Phys.  Plasmas  4, 1117  (1997)].  So,  tailoring  density  profiles  instead  of  using  flat  densities  is 
demonstrated  to  reduce  the  "seed"  mass  perturbation  amplitude  at  the  onset  of  the  exponential  RT 
growth. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF 

RESPONSIBLE  PERSON 

a.  REPORT 

unclassified 

b.  ABSTRACT 

unclassified 

c.  THIS  PAGE 

unclassified 

Same  as 
Report  (SAR) 

30 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


I.  INTRODUCTION. 

A  critical  issue  of  direct-drive  laser  fusion  program  is  mitigating  the  distortion  of  the 
imploded  pellet  due  to  the  Rayleigh-Taylor  (RT)  instability.1  The  common  approach  is  to  deal 
directly  with  the  initial  sources  of  the  instability.  Advances  in  improving  optical  smoothing  of  laser 
beams  and  developing  the  pellet  manufacturing  technologies  reduce  the  seed  for  the  ablative  RT 
instability  to  a  very  low  level.  Still,  even  a  small  perturbation  seed  could  be  large  enough  to 
preclude  the  achievement  of  high  gain  in  the  implosion. 

To  complement  the  two  above  direct  means  to  the  RT  seed  reduction,  consider  the  following 
approach.  Suppose  that  the  sources  of  mass  perturbations,  that  is,  the  levels  of  the  beam  and  target 
non-uniformity,  are  given.  Let  us  try  to  determine  the  structure  of  the  pellet  that  minimizes  the  seed 
for  the  RT  instability  produced  by  these  two  sources.  Here,  we  limit  ourselves  to  the  relatively 
simple  case  of  structuring  the  target  via  radial  density  tailoring.  The  high  gain  implosion  regime 
imposes  certain  requirements  on  the  thickness  and  the  areal  mass  of  the  pellet,  and  the  laser  pulse 
shape.  We  have  to  find  a  way  to  distribute  the  given  mass  within  a  layer  of  a  given  thickness  (that 
is,  the  way  the  density  profile  p(r)  should  be  tailored)  in  order  to  minimize  the  early-time  mass 
perturbation  amplitude.  In  general,  the  search  for  the  most  perturbation-resistant  pellet  structure 
could  include  tailoring  of  both  density  and  material  composition,  as  well  as  adjusting  the  laser  pulse 
shape  to  the  density/material  structure. 

It  is  known"'  that  the  fastest  classical  RT  mode  is  localized  near  the  accelerated  surface,  and 
its  growth  rate  is  not  sensitive  to  the  density  profile  in  the  fluid.  There  are  some  indications4  that 
density  tailoring  can  help  mitigate  the  target  distortion  due  to  the  exponential  RT  instability 
development,  but  this  is  not  the  effect  studied  here.  In  the  present  paper,  we  concentrate  on  the 
early-time  perturbation  growth  preceding  the  exponentiation  phase.  As  shown  from  recent  studies,5’ 
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the  evolution  of  perturbations  at  early  times  resembles  the  development  of  vorticity-driven 
Richtmyer- Meshkov  (RM)  instability  rather  than  the  acceleration-driven  RT  instability.  The  mass 
ablation  turns  out  to  be  a  sufficiently  strong  stabilizing  effect  to  stop  the  RM-like  perturbation 
growth  at  all  wavelengths  (in  contrast  with  the  RT  growth,  which  is  only  stabilized  at  sufficiently 
short  wavelengths  ).  Due  to  the  “rocket  effect”  caused  by  mass  ablation,  ’  mass  perturbations 
instead  of  growing  linearly  exhibit  decaying  oscillations.  These  oscillations,  discussed  in 
connection  with  ablative  RT  instability  ’  ,  were  first  observed  in  numerical  simulations  of  early-time 
mass  perturbation  growth.1’5  6  They  have  been  described  by  the  analytical  model  of  Ref.  7.  The 
initial  surface  roughness  and  laser  beam  non-uniformity  evolve  during  the  shock  transit  time  into 
mass  perturbations.  These  mass  perturbation  amplitudes  seed  the  RT  growth,  once  the  rarefaction 
wave  from  the  rear  side  of  the  target  reaches  the  ablation  front  and  the  subsonic  acceleration  of  the 
target  begins. 

The  role  of  density  tailoring  stems  from  the  possibility  to  invert  the  acceleration  of  the 
unstable  interface  during  the  shock  transit  time.  If  a  shock  wave  is  driven  into  increasing  density, 
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and  slows  down  as  it  propagates,  the  ablation  front  that  follows  the  shock  front  would  be 
decelerated,  too.  Hence,  local  gravity  is  in  the  same  direction  as  the  density  gradient,  and  small 
perturbations  tend  to  oscillate  (propagate  as  running  gravitational  waves  along  the  accelerated 
interface)  rather  than  growing  exponentially.  This  mechanism  contributes  to  the  oscillations 
additionally  to  the  “rocket  effect”  mentioned  above.  Consequently,  it  affects  the  early-time  evolution 
of  perturbations  by  increasing  characteristic  frequencies  of  their  oscillations.  For  a  given  level  of 
non-uniformity  supplying  the  initial  perturbations  (e.  g.,  lateral  variation  of  laser  intensity), 
increased  oscillation  frequencies  imply  lower  mass  perturbation  amplitudes  at  the  onset  of  the 
exponential  RT  growth. 

Use  of  tailored  density  profiles  for  stabilization  was  recently  proposed  to  mitigate  the  RT 
instability  development  in  imploding  Z-pinches. 10,11  As  far  as  the  perturbation  growth  is  concerned, 
the  main  difference  between  the  laser-driven  ablative  acceleration  and  the  electromagnetic  plasma 
acceleration  in  a  Z-pinch10  is  the  purpose,  and,  by  inference,  the  regime  of  acceleration.  In  a  Z- 
pinch  plasma  radiation  source  (PRS),  the  plasma  is  accelerated  in  order  to  convert  its  kinetic  energy 
into  thermal  at  stagnation,  and  finally  into  radiation.  Thus,  the  shock  transit  time  might  include  all  or 
most  of  the  implosion,  as  is  the  case  for  uniform  fill  or  thick  annular  cylindrical  loads.  A  shock- 
driven  acceleration  is  not  feasible  in  direct  drive  laser  fusion  -  too  much  pellet  preheat  prevents  fuel 
compression.  The  opportunity  to  affect  the  perturbation  growth  using  the  density  tailoring  is 
therefore  limited  to  the  early  time,  while  the  shock  and  compressibility  effects  are  still  relevant. 
Therefore,  the  ablation  drive  case  requires  a  much  more  detailed  study  of  the  perturbation  evolution. 

As  demonstrated  below,  the  radial  density  tailoring  affects  perturbations  whose  wavelengths 
are  of  order  of,  or  less  than  the  shell  thickness.  At  early  time,  while  the  aspect  ratio  is  large,  the 
spherical  pellet  is  nearly  planar  for  these  perturbations.  Therefore,  in  the  present  paper  we  study  the 
case  of  a  planar  geometry.  This  choice  somewhat  simplifies  both  analysis  and  numerics,  but  in  no 
way  reduces  the  generality  of  the  result.  The  stabilizing  effect  would  work  in  the  same  way  for 
spherical  geometry  as  it  does  for  planar  or  cylindrical10  cases. 

The  paper  is  structured  as  follows.  In  Section  II,  we  describe  the  dynamics  of  a  shock  wave 
driven  by  an  increasing  ablation  pressure  into  a  tailored  density  profile,  compare  the  1-D  simulation 
results  to  the  exact  self-similar  solution,  and  make  analytical  estimates  for  the  evolution  of 
perturbations.  Section  III  presents  a  numerical  study  of  the  perturbation  evolution  for  the  various 
types  of  density  tailoring,  and  Section  IV  concludes  with  a  discussion. 

II.  THEORY 

A.  Self-similar  unperturbed  flows 

Let  us  describe  one-dimensional  (1-D)  gas  dynamics  of  a  shock-piston  flow  generated  when 
a  shock  wave  is  driven  into  a  tailored  density  profile  by  a  pressure  that  either  remains  constant  or 
increases  with  time.  The  piston  could  be  visualized  as  an  ablation  front  for  laser  acceleration,  or  as 
the  magnetic  field/plasma  interface,  in  a  Z-pinch.10  We  will  find  how  the  piston  acceleration  (in  the 
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x  direction)  depends  on  the  relation  between  the  rate  of  pressure  increase  and  the  shape  of  the 
tailored  density  profile.  Typically,  the  plasma  represents  a  “heavy  fluid”,  accelerated  by  the  pressure 
of  a  “light  fluid”,  whether  it  be  a  low-density  ablated  plasma  or  a  massless  magnetic  field,  which 
gives  rise  to  the  RT  instability.  However,  the  presence  of  a  tailored  density  profile  in  a  shock-piston 
situation  could  reverse  the  direction  of  acceleration.  In  particular,  if  the  pre-shock  density  increases 
with  time  faster  than  the  post-shock  pressure,  the  interface  accelerates  opposite  to  the  propagation 
direction  of  the  shock  wave.  During  this  time,  the  propagating  shock  wave  keeps  bringing  new 
layers  of  the  initially  resting  gas  into  motion  in  the  positive  x  direction.  This  has  been  previously 
shown  heuristically,  using  a  0-D  snowplow  model.10  In  this  Section,  we  will  demonstrate  the  same 
in  a  more  rigorous  way,  with  the  aid  of  an  exact  self-similar  solution,  generalizing  that  obtained  in 
Ref.  12  and,  in  a  slightly  different  form,  in  Ref.  13.  Evolution  of  small  perturbations  on  top  of  this 
self-similar  flow  will  be  studied  below  -  analytically  in  Section  IIB  and  numerically  in  Section  III. 

Consider  a  half-space  x>0  ,  ,  ri  .  ,  , 

filled  with  an  ideal  gas  whose  density  profile  is  tailored 

according  to  a  power  law:  p(x)  =  pQ(x  /  A)x  ,  where  Po  and  A  are  dimensional  normalization 

constants,  and  the  constant  dimensionless  exponent  X  >  ~1  determines  the  shape  of  the  profile.  Let 

a  piston,  which  moves  according  to  the  power  law 

xp(t)  —  L{t)  =  g0tm  t  m  (1) 

(here  again,  g0  is  a  dimensional  normalization  constant,  m  is  a  dimensionless  exponent)  drive  a 


shock  wave  into  the  gas.  In  the  reference  frame  where  the  piston  is  at  rest,  the  1-D  motion  of  the 


shocked  fluid  is  described  by  the  equations  of  continuity,  adiabaticity,  and  motion: 
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where  u,  p,  p,  and  yare  the  x- velocity,  density,  pressure,  and  adiabatic  exponent  of  the  shocked 


gas,  respectively.  For  m  >  1,  both  the  shock  and  the  piston  accelerate  with  time,  so  that  local 
gravity  has  negative  direction,  from  the  accelerated  fluid  to  the  piston.  For  m  <  1,  the  shock  and  the 
piston  decelerate,  and  the  local  gravity  is  positive,  directed  from  the  piston  to  the  fluid. 

If  the  pre-shocked  gas  is  cold,  so  that  the  shocks  wave  is  strong  at  all  time,  then  Eqs.  (2)- 
(4)  admit  an  exact  self-similar  solution.  Introducing  a  self-similar  variable,  V  =  x  / L(t),  we  use  the 


self-similarity  ansatz 

u(x,t)  =  Ut)U(Jt)  =  gQt"l~1U  (77), 

p(x,t)  =  p0L(t)x  N  (jj)  =  p0(g0/m)xtmxN(ri),  (5) 

p(x,t)=  p0L(t)x  L2P(rj)=  p0(gx+2/rnx)t2m-2+mxP(j1), 
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where  U (rj),  N(rj), and  P{tj)  .  . 

are  dimensionless  functions  representing  the  self-similar  profiles  of 

velocity,  density,  and  pressure,  respectively.  With  the  aid  of  (5),  the  gasdynamic  equations  (2)-(4) 

are  reduced  to  a  system  of  two  ordinary  differential  equations,  which  has  to  be  solved  numerically, 

as  described  in  Appendix  I. 

As  seen  from  Eq.  (5),  this  solution  is  generated  if  the  pressure  at  the  piston  is  related  to 
either  distance  traveled  by  the  piston  or  to  time  by  a  power  law,  e.  g.,  pp{t)  <=<  f .  Comparing  this 


to  (5),  we  find  the  value  of  m  that  corresponds  to  <7  for  a  given  %: 
(7  +  2 


m  = 


Z  +  2 


(6) 


Equation  (6)  illustrates  that  the  shock-piston  flow  accelerates  the  fluid  when  the  pressure  at  the 
piston  rises  faster  than  the  density  in  the  unperturbed  fluid:  m  >  1  when  o  >  %.  Otherwise,  when 
<7  <  x  ,  the  fluid  is  decelerated  (m  <  1).  For  instance,  if  the  pressure  rises  linearly,  a  =  1 ,  then  the 
plasma  layer  of  parabolic  density  (%=  2)  will  be  decelerated:  m  =  3/4.  If  the  same  linearly  rising 
pressure  drives  a  uniform  density  profile  (%=  0),  then  the  pre-shock  density  increases  slower  than 
the  post-shock  pressure  does,  which  means  acceleration:  m  =  3/2. 

The  self-similar  profiles  of  post-shock  density  and  temperature  in  an  ideal  gas  with  y  =  5/3 
for  these  two  cases  are  shown  in  Fig.  1.  The  density  profile  near  the  piston,  at  77  «  1,  has  the 
power-law  asymptotic  shape 

%m(y- 1)  -  2  (m—  1) 


A (77)  ~  77s ,  where  s  = 


(7) 


m(y  +x  +2)  -2 

which  is  typically  sharper  than  the  pre-shock  profile.  Since  the  pressure  at  V  n  ^  must  be  of  order 
unity  (PQl)  ~  r7° ),  we  find  for  the  self-similar  temperature  profile:  ©(77)  =  P(ij)  /  N (77)  ~  77  ' .  For 
the  above  two  cases,  c  = 1 ,  X  =  0  and  2,  we  find  from  (7):  s  =  -2/7  and 2/3  ,  respectively.  Then 
the  temperature  vanishes  and  the  density  diverges  near  the  accelerated  piston  (x  =  0,  m  =  3/2),  and 
vice  versa,  the  temperature  diverges  and  the  density  vanishes  near  the  decelerated  piston  (X=2,  m 
=  3/4),  cf.  Figs  .1  (a),  (b)  This  is  not  surprising.  The  entropy  of  the  fluid  particles  adjacent  to  the 
piston  is  determined  at  1  n  when  the  shock  wave  is  formed.  For  m  =  3/2,  the  shock  velocity  in 


this  limit  tends  to  zero,  and  so  do  the  post-shock  temperature  and  pressure.  The  isentropic 
compression  that  follows  can  make  the  pressure  in  these  fluid  particles  finite  only  via  infinite 
increase  of  density.  Conversely,  for  m  =  3/4,  the  shock  velocity  in  the  limit  t  n  0  diverges,  and  so 
does  the  post-shock  temperature.  Due  to  the  subsequent  isentropic  compression,  the  temperature 
there  remains  infinite  afterwards. 

Figure  2  presents  a  comparison  of  the  self- similar  density  profile  calculated  for  the  case  of 
decelerated  piston  ( o  =  1 ,  X  =  2,  m  =  3/4)  with  those  obtained  in  our  simulations  for  similar 
conditions.  The  simulated  density  profiles  are  seen  to  remain  self-similar  during  the  shock  transit 
time.  The  density  gradient  near  the  ablative  piston  in  simulation  is  less  steep  than  in  the  analytical 
solution,  because  the  temperature  rises  even  faster  than  in  Fig.  1(a),  the  ablation  pressure  being 
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maintained  primarily  through  the  conductive  heat  flow  from  the  corona.  Also,  in  the  computations 
the  density  at  x  =  0  is  finite.  Nevertheless,  the  observed  self- similarity  supports  our  prediction  that 
the  ablation  front  in  this  case  is  decelerated,  as  shown  below  in  Section  III. 

B.  Simple  linear  stability  analysis 

Now  we  assume  that  the  piston  driving  the  self-similar  flow  is  an  ablation  front.  Then  the 
driving  pressure  is  created  by  a  mass  ablation  through  this  front,  so  that  we  deal  with  a  “leaky 
piston”  rather  than  a  solid  piston  or  a  contact  interface.  Still,  the  self-similar  solution  of  Section  II.A 
could  be  valid,  once  the  ablation  front  is  identified  with  some  surface  77  =  rja  =  const  rather  than 
with  the  surface  77  =  0.  Then  the  ablation  velocity  va  ,  like  all  the  characteristic  velocities,  is 


proportional  to  L  =  mL/t ,  where  L  is  the  distance  traveled  by  the  piston,  see  Eqs.  (1),  (5).  We 
introduce  a  dimensionless  constant  v  characteristic  of  the  ablation  rate: 

V"  (8) 


mL/t 


Let  us  describe  perturbations  of  the  ablation  front  using  a  discontinuity  model  identical  to 
that  of  Piriz  et  al.9  The  model  involves  another  dimensionless  characteristic  of  the  ablation  front, 
the  density  ratio  rD  =  p2  /  pt ,  where  the  subscripts  1  and  2  refer  to  the  dense  shocked  plasma  ahead 

of  and  to  the  low-density  plasma  behind  the  ablation  front,  respectively.  This  parameter  is  also 
supposed  to  be  constant  in  time.  Typically,  rD«  1 ,  so  that  the  Atwood  number  across  the  ablation 
front,  A  =  (1  -  rD )  /(I  +  rD) ,  is  close  to  unity.  We  modify  the  boundary  conditions  given  in  Ref.  9 
for  the  case  of  time-dependent  acceleration  and  ablation  velocity,  and  require  all  the  perturbations  to 
be  localized  near  the  ablation  front  or  propagating  from  it,  upstream  or  downstream.  Then  the 
solution  of  the  perturbation  problem  contains  no  perturbations  incident  on  the  front,  only  waves 
outgoing  from  it.  As  shown  in  Appendix  II,  the  boundary  conditions  in  this  case  are  reduced  to 

— 8x  +  kv  8x  -  Sv  =  0,  (9) 

7  .  ci  a  ci  x  7 

at 


—8vx  +  - — —kva8vx  +A  - — —lev:  - gk 
at  1+  rD  |_  rD 


8xa  =0, 


(10) 


where  8xa  is  the  displacement  amplitude  of  the  ablation  front,  8v  v  is  the  perturbation  of  the  axial 


velocity  at  the  ablation  front,  k  is  the  perturbation  wave  number.  The  boundary  condition  (9)  at  the 
ablation  front  is  well  known  in  the  literature,14  the  boundary  condition  (10)  is  derived  in  Appendix 
II  from  the  results  of  Ref.  9  under  the  above  assumptions. 

Introduce  a  new  dimensionless  independent  variable  [not  to  be  confused  with  one  of  the 
dependent  variables,  self-similar  temperature  ©(77)] 

6  =  kLit)  =  kg  Qtm  /  m  .  (11) 

Substitution  of  (8),  (11)  and  kg(t )  =  m(m- 1)0  It 2  into  Eqs.  (9),  (10)  reduces  this  system  to  a 
single  equation  for  8xa : 
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8x  =  0 . 


Solutions  of  this  equation  are  expressed  via  the  confluent  hypergeometric  functions.  In  particular, 
the  solution  regular  at  t  =  0  and  satisfying  the  initial  condition  8xa  (t  =  0)  =  8x0  is  given  by 

gk  =  eJ-vf_g_+ -gdj  1  (13) 


-j— '  =  exP  -V  - - +  -=  0  - - -  1-- - ,(13) 

&o  [  l1+r°  \r°  J  J  {  m  L  v  J  m  \lr»  J 

where  ^,(<3,0,2)  is  the  Kummer’s  function.15 

To  analyze  the  properties  of  this  solution,  let  us  start  with  the  simple  case  of  negligible  effect 
of  mass  ablation,  obtained  from  (13)  in  the  limit  v  — >  0 .  In  this  limit,  the  Kummer’s  function  is 
reduced  to  the  Bessel  function  (Ref.  15,  Eq.  13.5.13): 


m—  1  )  (1  -  in)  AG 


(1  -m)A0 


Since  A  >  0,  the  evolution  of  perturbations  is  determined  primarily  by  the  sign  of  1-  m,  that  is,  by 
the  direction  of  acceleration  of  the  interface.  If  m  <  1,  so  that  the  interface  is  decelerated,  and  there 
is  no  reason  for  the  RT  instability  to  develop,  then  the  argument  of  the  Bessel  function  in  (14)  is 
real.  At  late  time  ( 0  » 1) ,  the  perturbations  oscillate: 


cos  2 


8x  1  [  m- 1  ]  (1  -m)A0  4m  [_  1(1  -m)A6  n  n\ 

8x 0  yjn  ^  m  J\_  m  \  V  m  2 m  4 

The  amplitude  of  oscillations,  as  seen  from  Eq.  (15),  is  slowly  growing.  A  qualitative  explanation 
of  this  is  found  in  Ref.  10:  amplitude  of  a  pendulum  oscillating  in  a  slowly  decreasing  gravitational 
field  increases. 

For  m  >  1,  the  interface  is  accelerated,  and  therefore,  is  RT  unstable.  In  this  case,  in  Eq. 
(14)  we  have  a  modified  Bessel  function,  which  describes  exponential  growth  instead  of 
oscillations  at  late  time: 


2 -JtF 


m  —  1  I  (m  - 1)  A  0 


exp  2 


l(m-\)A0 


For  the  particular  case  of  constant  acceleration  (m  =  2),  Eq.  (14)  yields 

coshQ  Agkt),  (17) 

as  it  should  be.  The  perturbation  growth  described  by  Eq.  (14)  for  the  cases  of  m  =  3/4  and  m  =  3/2 
is  shown  in  Fig.  3(a). 

Introduction  of  mass  ablation  (finite  v)  changes  the  picture  qualitatively.  Now  the 
asymptotic  behavior  of  solution  (13)  at  late  time  is15 

^Xa  _  .  Q~T>0  _,,f  2  ■  A  1  ns) 


—2-  ~  const  ■  0  2m  Rel  exp  -  v  - +i  —=  \6  +  i 
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The  solution  describes  decaying  oscillations  for  any  sign  of  acceleration.  It  is  qualitatively  similar  to 
the  particularly  simple  case  of  zero  acceleration  [  m  =1,6  =  kv0t,  where  v0  for  this  case  is  an 
appropriate  replacement  of  g0  in  the  definition  (1)]: 


K 

8x„ 


exp 


2  vkvQt 
1  +  rn 


cos 


f  vAkvQt  ^ 

jo 


(19) 


cf.  Ref.  7.  The  evolution  of  perturbations  described  by  Eq.  (13)  for  the  cases  of  m  =  3/4  and  m  = 
3/2  and  v=  0.05  is  shown  in  Fig.  3(b).  The  only  difference  between  the  cases  of  acceleration  (m  = 
3/2)  and  deceleration  ( m  =  3/4)  is  that  in  the  former  case  the  perturbation  amplitude  reaches  a  higher 
peak  value  and  oscillates  at  a  lower  frequency. 

Stabilization  of  all  modes  in  all  acceleration  regimes  is  explained  with  the  aid  of  the 
dispersion  relation  obtained  in  Ref.  9  (self-consistent  derivations  of  similar  relations  are  found  in 
Refs.  8,  16).  The  exponential  growth  rate  for  ablative  RT  instability  is  found  to  be  equal  to 


r  = 


\l/2 


2.,  2 


Agk  —  ——k  v 


2kVa 

1  +rn 


(20) 


This  expression  follows  from  Eqs.  (9),  (10)  if  the  ablation  velocity  and  acceleration,  va  and  g,  are 
assumed  constant.  For  the  wavenumbers  exceeding  a  cut-off  value  kc ,  the  instability  is  fully 


suppressed.  Equation  (20)  demonstrates  that  the  cut-off  wavelength  is  determined  by  the  balance  of 
the  two  terms  in  parentheses,  k  va  /  rD  >  gk  ,  rather  than  by  the  balance  of  the  whole  term  in 

parentheses  and  the  second  terms  in  the  right-hand  side  of  Eq.  (20),  kva  >  gk  ,  since  rD  « 1.  For 
k  exceeding  the  cut-off  value,  Eq.  (20)  predicts  exponentially  decaying  oscillations.  Qualitatively, 
this  stabilization  mechanism  is  described  in  Ref.  9.  It  is  based  on  the  fact  that  the  ablation  front  is  an 
isotherm.  When  the  front  moves  into  a  higher  temperature  area,  its  temperature  therefore  does  not 
increase,  but  the  temperature  gradient  near  it  does.  Increased  temperature  gradient  speeds  up  the 
mass  ablation,  producing  the  “rocket  effect”  responsible  for  the  oscillations  that  tend  to  decrease 
the  perturbation. 

For  our  self-similar  acceleration  regimes,  k2v2  t2'"~2 ,  whereas  gk  °c  tm~2 .  Since  m  is 
positive,  k2v2  grows  with  time  faster  than  gk.  If  g  is  positive  ( m  >  1,  acceleration),  then  eventually 
k  2v]  /  rD  exceeds  gk,  and  the  exponential  perturbation  growth  evolves  into  decaying  oscillations. 


Otherwise,  if  g  is  negative  ( m  <  1,  deceleration),  we  have  decaying  oscillations  from  the  start,  see 
Eq.  (20).  This  is  exactly  what  is  described  by  Eqs.  (13),  (18),  (19),  and  illustrated  by  Fig.  3(b). 

To  summarize,  we  have  predicted  decaying  oscillations  instead  of  exponential  perturbation 
growth  for  all  wavelengths  and  all  acceleration  regimes  subject  to  self-similarity  restriction.  This  is  a 
physically  meaningful  conclusion  that  could  be  tested  in  simulation  and  experiment.  In  contrast  with 
that,  neither  the  prediction  of  the  exponential  decay  of  perturbation,  nor  the  calculated  rates  of  decay 
should  be  taken  very  seriously.  The  reason  is  that  our  simplifying  assumption  that  no  external 
perturbations  come  to  the  ablation  front  as  incident  waves  [see  Eq.  (A17)]  is  not  fully  justified. 


8 


Actually,  the  plasma  that  reaches  the  ablation  front  has  first  passed  through  the  rippled  shock  front. 
The  plasma  particles  therefore  bring  to  the  ablation  front  the  perturbations  of  vorticity  left  after  the 
shock  passage.  The  incident  perturbations  also  decay  with  time,  but  much  slower  than 
exponentially.  For  instance,  in  the  case  of  zero  acceleration  (m  =  1),  the  perturbations  oscillate  and 
decay  as  Bessel  functions  of  time,  that  is,  as  t~112  or  t~312 ,  depending  on  whether  the  shock  wave  is 
strong  or  not.  If  the  ablation  front  is  RT  unstable,  then  this  power-law  source  of  external 
perturbations  does  not  change  the  exponential  growth  rate.  This  is  why  the  incident  perturbations 
could  be  safely  neglected  when  the  expressions  for  the  growth  rate  like  (20)  are  derived  (if  the 
resulting  growth  rate  is  positive).  If  there  is  no  exponential  RT  growth,  this  approximation  is  no 
longer  justified.  Therefore,  the  exponential  decay  rates  estimated  above  have  little  physical 
significance.  Still,  the  corresponding  oscillation  frequencies  are  physically  meaningful  because  they 
describe  the  immediate  pressure  response  to  the  perturbation  of  the  ablation  front  due  to  the  “rocket 
effect'’.9 

The  complicated  early-time  evolution  of  perturbations  generates  the  seeds  for  the  exponential 
RT  instability  development  that  starts  after  the  shock  transit  time.  We  are  interested  in  designing  our 
target  so  that  the  given  external  perturbations  (e.  g.,  lateral  non-uniformity  of  the  laser  beam) 
generate  as  small  mass  perturbation  seeds  as  possible.  This  is  where  the  tailored  density  profiles  can 
help.  Tailoring  the  initial  density  profile  makes  g  negative  from  the  start.  As  shown  above  [e.  g., 
see  Fig.  3(b)],  this  increases  the  frequency  of  oscillations  of  the  ablation  front  and  decreases  the 
peak  perturbation  amplitude.  We  expect  this  effect  to  help  decrease  the  perturbation  amplitudes  at 
the  time  when  the  RT  instability  starts  to  develop.  Simulation  results  presented  in  the  next  Section 
demonstrate  that  this  indeed  could  be  done. 

III.  NUMERICAL  SIMULATION 

The  simulations  were  performed  in  two  dimensions  (2-D)  using  the  FAST  hydrocode.19  This 

20 

code  has  proved  to  be  very  robust  over  a  wide  range  of  problems,  including  reactive  flows, 

21  22 

turbulent  jets,-  flows  about  ship  superstructure,"-  and  RT  instabilities  in  inertial  confinement  fusion 

23 

targets."  In  the  studies  of  the  early-time  perturbation  growth  in  laser  targets,  the  FAST  code  has 
proved  effective  at  reproducing  experimental  observations  of  laser  imprinting-  and  at  predicting  the 
oscillations  due  to  the  “rocket  effect”.1'6  The  code  includes  flux-corrected  transport  (FCT) 
advection,-  inverse  bremsstrahlung  laser  energy  absorption  (multiple  ray  trace),  and  Spitzer-Harm 
thermal  conduction.26  Radiation  transport  and  tabulated  equation  of  state  (EOS)  are  also  included  in 
the  code.  These,  however,  have  not  been  invoked  in  the  simulations  reported  below.  Our  goal  is  to 
demonstrate  the  hydrodynamic  effect  of  density  tailoring  alone,  without  any  contribution  of 
smoothing  due  to  radiation  transport  or  EOS  effects.  We  simulate  ablative  acceleration  of  CH 
(plastic)  targets. 

We  compare  four  target  structures  that  differ  mainly  by  their  density  profiles,  see  Fig.  4(a). 
[Note  that  in  the  simulations,  the  laser  radiation  comes  from  the  positive  x  direction,  accelerating  the 
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target  to  negative  x  (the  coordinate  x  emulates  r  in  spherical  geometry),  a  different  sign  from  that 
chosen  in  Section  II.  This  should  not  cause  any  confusion,  since  all  the  results  of  Section  II  are 
invariant  under  the  substitution  x  <->  —x .  ] 

All  the  targets  have  the  same  areal  mass  68.2  pmg/cm  .  Targets  1  to  3  have  also  the  same  total 
thickness,  140  (im.  The  front  part  of  target  1  approximates  the  parabolic  density  profile  of  Section 
II  (x  =  2).  Density  profile  starting  from  p  =  0  is  not  realistic  in  target  design,  so  that  the  actual 
density  profile  of  target  1  is 

P(x)  =  Po  +  (Pi-Po||-  ,  (21) 

where  p0  =  0.05  g/cm3  (a  foam  of  this  low  density  could  be  manufactured;  this  value  is  just  about 
the  critical  density  for  KrF  laser  radiation),  p,  =1.07  g/cm3  (solid  plastic  density),  A  =  120  pm, 
followed  by  a  20  pm  payload  of  solid  plastic  [we  also  did  some  runs  for  the  initial  density  profile 
(21)  without  adding  the  payload].  In  target  2,  the  areal  mass  contained  in  the  profile  (21)  is 
uniformly  distributed  over  the  same  thickness,  120  pm,  which  corresponds  to  a  uniform  density 
p  =  0.39  g/cm 3 ,  and  followed  by  the  same  20  pm  payload.  In  target  3,  the  total  mass  68.2 
pm^g/cm3  is  uniformly  distributed  with  the  constant  density  p  =  0.487  g/cm 3 ,  over  the  same  total 
thickness,  120  pm  +  20  pm  =  140  pm.  Finally,  the  solid  plastic  target  4  has  the  same  total  mass, 
and,  therefore,  is  63.7  pm  thick,  that  is,  thinner  than  the  other  three. 

The  targets  are  driven  by  a  0.248  pm  KrF  laser  radiation,  as  in  Nike  laser  at  the  Naval 
Research  Laboratory."  Two  shapes  of  laser  pulse  are  used.  One  is  selected  to  provide  in  target  1  a 
linear  increase  of  pressure  from  0  to  4  Mbar  during  the  shock  transit  through  the  density  profile 
(21),  that  is,  in  the  first  4  ns  [Fig.  4(b)],  after  which  the  laser  intensity  is  held  constant.  This  pulse 
shape,  as  expected  from  the  results  of  Section  II  and  demonstrated  in  Fig.  5  below,  at  early  time  (< 
4  ns)  produces  a  decelerating  ablation  front  in  target  1,  and  an  accelerated  ablation  front  in  targets  2 
to  4.  The  other  laser  pulse  shape  studied  here  is  a  constant  incident  laser  intensity,  which 
corresponds  to  the  ablative  pressure  slowly  decreasing  from  2.75  Mbar  to  1.75  Mbar  at  10  ns  due 
to  increased  absorption  in  the  plasma  corona.  Here  again,  the  ablation  front  in  target  1  at  early  time 
must  be  decelerating.  In  other  targets,  however,  we  expect  zero  acceleration,  or  even  a  small 
deceleration  (due  to  decreasing  pressure).  Of  course,  the  theory  of  Section  II  does  not  extend  to  the 
situation  after  the  shock  wave  breaks  out  at  the  rear  surface  of  targets  1,  3,  4,  or  is  reflected  from 
the  payload  boundary  of  target  2.  There  is  no  analytical  theory  available  to  describe  perturbation 
growth  during  the  reverberation  of  shock  and  expansion  waves,  which  eventually  leads  to  the 
subsonic  acceleration  regime.  Here,  we  have  to  rely  on  the  simulation  results. 

To  facilitate  comparison  between  various  cases,  we  introduce  a  constant  single-mode  lateral 
variation  of  the  incident  laser  intensity.  Evolution  of  perturbations  studied  below,  therefore, 
corresponds  to  laser  imprinting  of  mass  perturbations.  As  explained  in  Ref.  6,  the  early-time 
perturbation  growth  due  to  the  surface  roughness  is  physically  equivalent  to  that  produced  by  the 
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laser  imprint.  The  difference  is  only  in  the  initial  conditions:  in  the  former  case,  we  create  initial 
amplitude  of  the  perturbation,  whereas  in  the  latter  we  supply  the  perturbation  with  initial 
momentum.  The  perturbation  growth  that  follows  is  governed  by  the  same  physics  and  could  be 
scaled  from  one  case  to  another.  With  tailored  density  profiles,  the  same  amplitude  of  surface 
ripples  would  imply  a  20-fold  difference  in  the  initial  mass  variation  between  targets  1  and  4. 

Below,  we  have  all  the  targets  initially  uniform,  driven  by  laser  pulses  with  the  same  lateral 
intensity  variation,  so  that  the  comparison  between  the  targets  is  deemed  fair.  For  each  target/laser 
pulse  shape  combination,  we  run  the  code  with  single-mode  perturbation  wavelengths  of  15  (im,  30 
pm,  and  60  pm. 

Figure  5  is  plotted  for  targets  1  (a)  and  2  (b)  driven  by  a  linearly  increasing  ablative 
pressure.  It  shows  the  x-t  trajectories  of  fluid  particles,  that  give  an  idea  about  the  motion  of  the 
ablation  front.  We  see  that  at  early  time  the  ablation  front  in  target  1  is  indeed  decelerated,  as 
predicted  in  Section  II.A.  Similarly,  we  demonstrate  that  for  all  the  other  targets,  the  ablation  front 
is  accelerated  by  the  linear  pressure  increase  from  the  beginning  of  the  laser  pulse. 

Figures  6  to  8  compare  the  evolution  of  mass  perturbations  between  targets  1  to  4,  all  of 
them  driven  by  a  laser  pulse  with  linearly  increasing  ablative  pressure,  and  a  0.2%  rms  single-mode 
lateral  intensity  variation,  X  =  15  pm,  30  pm,  and  60  pm,  respectively.  Here,  the  early  time  (shock 
transit)  corresponds  to  the  first  4  ns,  whereas  the  “late  time”,  when  the  RT  instability  of  acceleration 
develops  exponentially,  starts  somewhere  between  6  and  8  ns. 

For  all  cases,  the  qualitative  behavior  agrees  with  the  predictions  of  Section  II.B.  The  mass 
perturbations  tend  to  oscillate,  no  matter  whether  the  ablation  front  is  accelerated  or  decelerated, 
although  for  the  latter  case  (target  1)  the  frequency  of  oscillations  is  notably  higher,  and  the 
amplitude  lower.  We  even  observe  some  [albeit  non-exponential,  cf.  Fig.  3(b)]  decrease  in  the  mass 
variation  amplitude  for  the  shortest  perturbation  wavelength  X  =  15  pm,  see  Fig.  6.  Although  we 

have  not  studied  the  scaling  of  oscillation  frequencies  with  the  perturbation  wavelength  [in  terms  of 
the  discontinuity  model  of  Ref.  9  and  Appendix  II,  even  the  scaling  of  rD  with  X  is  not  known  for 

tailored  density  profiles],  one  can  naturally  expect  longer  perturbation  wavelengths  to  correspond  to 
lower  frequencies.  This  is  indeed  the  case.  During  the  shock  transit  time,  we  have  multiple 
oscillations  in  target  1  for  X  =  15  pm,  one  full  cycle  for  X  =  30  pm,  and  half  cycle  for  X  =  60  pm. 
For  all  the  other  cases,  we  have  no  more  than  a  half  cycle. 

Eventually,  the  oscillatory  behavior  evolves  into  an  exponential  growth  for  all  targets.  The 
RT  exponential  growth  rates  T  =  *Jgk  arc  essentially  the  same  for  all  targets,  since  they  are  driven 
by  the  same  ablative  pressure,  have  the  same  mass,  and  therefore  experience  the  same  acceleration. 
There  is,  however,  a  significant  difference  between  the  targets  in  the  amplitudes  of  these 
exponentially  growing  modes.  This  is  precisely  the  difference  in  the  seeds  for  RT  discussed  above. 
Firstly,  we  notice  that  the  perturbations  in  the  solid  target  4  are  consistently  higher  than  in  the  other 
three.  This  clearly  illustrates  the  beneficial  effect  of  decreasing  the  average  density  of  the  target.  The 
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effect  has  nothing  to  do  with  radiative  smoothing  of  perturbations  (recall  that  to  emphasize  this 
point,  we  have  no  radiation  transport  here),  it  is  purely  hydrodynamic,  as  explained  in  Ref.  6. 
Secondly,  the  perturbations  in  target  1  containing  a  tailored  density  profile  are  consistently  lower 
than  in  the  other  three.  Evolution  of  perturbations  in  targets  2  and  3,  4  is  differently  affected  by  the 
shock/rarefaction  reverberations  (a  material  interface  in  target  2  causes  an  additional  shock 
reflection),  and  thereby  phase-sensitive:  target  2  produces,  a  higher  or  a  lower  RT  seed  than  target  3 
for  A  =  30  (im  and  A  =  60  pm,  respectively,  whereas  for  A  =  15  pm  both  perform  similarly.  In 
contrast  with  this,  target  1  always  produces  the  lowest  RT  seed,  which  agrees  with  the  results  of 
Section  II.B.  The  stabilizing  effect  of  radial  density  tailoring  is  generally  more  effective  for  shorter 
wavelengths  (at  least  a  half  cycle  of  oscillations  has  to  be  completed  during  the  shock  transit  time, 
cf.  Figs.  6-8),  but  it  does  not  otherwise  discriminate  between  the  perturbation  wavelengths.  In  other 
words,  the  effect  is  sufficiently  robust. 

This  conclusion  is  confirmed  by  inspection  of  Figs.  9  to  11,  plotted  for  a  constant  intensity 
laser  drive,  all  the  other  conditions  being  the  same  as  in  Figs.  6  to  8,  respectively.  Here, 
acceleration  of  targets  2  to  4  at  early  time  is  very  low.  We  observe  oscillations  due  to  sonic  waves 
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propagating  in  the  lateral  direction,  as  in  the  other  cases  of  the  RM-like  instability  development. 
Otherwise,  the  situation  is  quite  similar  to  that  presented  in  Figs.  6  to  8:  the  tailored  density  target  1 
is  consistently  the  best,  the  uniform  solid  target  4  is  the  worst,  whereas  the  double-slab  target  2  and 
the  low-density  uniform  target  3  are  somewhere  in  between,  their  relative  performance  depending 
on  the  perturbation  wavelength. 

Figures  12,  13  that  show  density  perturbation  maps  of  p  -M  ,  where  M  is  the 


average  of  the  density  over  y  at  a  given  x,  in  targets  1  and  3  at  various  instants  of  time,  help  to  see 
why  the  mass  perturbations  are  lower  in  a  target  with  a  tailored  density  profile.  In  target  1,  at  early 
time  the  density  perturbation  is  spread  over  a  wide  plasma  volume.  The  oscillations  tend  to 
distribute  the  extra  density  more  uniformly.  The  bubble-spike  structure,  distinct  for  a  developed 
single  mode  RT  instability,  is  only  formed  at  late  time.  In  target  3,  the  bubble-spike  structure, 
developing  from  the  start,  oscillates  as  a  whole.  During  the  oscillation  half  cycle,  the  bubbles  and 
the  spikes  exchange  places,  but  they  are  always  present.  Consequently,  the  mass  perturbations  are 
more  localized  in  target  3  than  in  target  1,  and  it  counts  when  the  RT  growth  (which  is  precisely  the 
growth  of  localized  mass  perturbations)  starts. 


IV.  CONCLUSIONS 

Radial  tailoring  of  initial  density  profiles  in  laser  targets  inverts  the  acceleration  of  the 
ablation  front  at  early  time,  thereby  increasing  the  frequency  and  decreasing  the  amplitude  of  its 
oscillations.  This  effect  reduces  the  seed  for  the  exponential  RT  instability  growth  generated  by  a 
given  lateral  non-uniformity  in  the  laser  drive  in  a  target  of  given  mass  and  thickness.  A  significant 
reduction  in  the  seed  mass  perturbation  amplitude  in  a  target  with  a  tailored  density  profile  in  our 
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simulations  was  found  for  all  perturbation  wavelengths,  compared  to  a  uniform  target  of  the  same 
mass  and  thickness,  at  the  level  of  at  least  5  or  more. 

The  stabilizing  effect  of  density  tailoring  is  purely  hydrodynamic.  In  a  way,  it  complements 
a  related,  but  different,  stabilizing  hydrodynamic  effect  of  decreasing  the  initial  density  of  the  target, 
which  also  reduces  the  seed  for  the  RT  instability.6  If  the  target  design  includes  low-density  outside 
foam  layers  and  more  dense  inner  plastic  or  DT  layers,  then,  as  demonstrated  in  our  simulations,  a 
smooth  transitions  from  low  to  high  initial  density  provides  a  significant  reduction  in  the  RT  seed 
compared  to  the  case  of  two  adjacent  layers  of  different  densities. 

Our  present  results  suggest  the  next  step.  A  tailored  density  profile  should  be  incorporated 
into  a  high-gain  target  design  for  direct-drive  laser  fusion,  and  an 

implosion  of  such  a  target  should  be  simulated  with  all  the  relevant  physics  included,  to  see  whether 
we  could  indeed  expect  a  robust  reduction  in  the  net  mass  variation  at  stagnation.  This  work  is  now 
in  progress. 


Acknowledgments 

We  are  grateful  to  Drs.  S.  E.  Bodner  and  J.  P.  Dahlburg  (NRL)  for  stimulating  discussions. 
One  of  us  (N.  M.)  wishes  to  thank  Dr.  Bodner  and  Laser  Plasma  Branch  at  NRL  for  their  support 
and  hospitality  during  the  performance  of  this  research.  The  work  was  sponsored  by  the  U.S. 
Department  of  Energy  through  the  contract  for  the  Naval  Research  Laboratory. 


13 


Appendix  I.  Seif-similar  solution  for  the  shock-piston  flow  with  a  tailored 
density  profile. 

Substituting  the  ansatz  (5)  into  Eqs.  (2)-(4)  and  introducing  a  new  self-similar  variable  by 
S(rf)  =  yP{ri)  /{[TJ-U  ( ti)]N(t) ) } ,  we  obtain  the  following  system  of  ordinary  differential  equations 


for  the  self- similar  variables: 


(u-m 

(U-rt) 

dU  _ 
dr\ 


d\nN  dU 

-+ - +X=0 , 


dri 

dhiS 

dri 


dri 

dU 


-+y— +1-  — =  0, 
di]  m 


(2m  -  2  +  mx)S  +  y(m  —  1)(C7  + 1) 
ym(U  +S-tj) 


(Al) 


(A2) 

(A3) 


Substituting  (A3)  into  (A2),  we  obtain 

dlnS  _  in(  J  +  1)3  +  [/( in  - 1)  +  2  -  m ]U  +  ( in  -  2)7]  +  y(in  -  1) 

~dri  _  m(U  +S-rj)(T]-U)  '  ( 

The  two  equations  (A3)  and  (A4)  must  be  solved  numerically.  Note  that  N{ri)  does  not  enter  Eqs. 
(A3)-(A4).  In  a  general  case,  once  a  solution  of  Eqs.  (A3),  (A4)  is  found,  N(rj)  could  be  expressed 
via  1/(77)  and  S(ri)  from  Eqs.  (A1)-(A2)  as 

",{X+l)  m(n-l)+2 

N  =  const  ■  S",(r+l)~2  ( n-U  )m  r+i>-2 

(A5) 


(the  adiabaticity  integral),  where  the  constant  is  determined  from  the  boundary  conditions  at  the 
shock  front.  In  a  particular  case  of  m  =2/(y+  1)  (for  y=  5/3,  this  is  our  case  of  decelerated  piston), 

the  integral  (A5)  cannot  be  used,  but  instead  Eq.  (A2)  is  immediately  integrated  to  obtain 
S  (r/)  =  const  ■  (77  -  U  fy .  Then  this  expression  for  S(  i])  is  substituted  into  Eq.  (A3),  which  must  be 

integrated  together  with  (Al).  In  either  case,  we  have  to  solve  two  coupled  ordinary  differential 
equations. 

The  boundary  condition  at  the  piston,  77  — >  0 ,  is  U  — >  0 ,  hence,  S  — >  00  .  Then  we  find 
from  (A3)  that  near  77  =  0 

. .  2m  — 2  + my  , .  „ 

U  ~  - - —77.  (A6) 

ym 

After  the  shock  formation,  the  gas  near  the  piston  is  compressed  (decompressed)  if  the  coefficient  in 
the  right-hand  side  of  Eq.  (A6)  is  negative  (positive).  Substituting  (A6)  into  (A4),  we  find  an 
asymptotic  expression  for  S(  i]): 

__W+I)_ 

S  ~SQ-ri  (2+r+x)m~2,  (A7) 

where  S0  is  a  constant  to  be  determined.  Substituting  (A6)  and  (A7)  into  the  integral  (A5),  we 


derive  Eq.  (7). 

The  boundary  conditions  at  the  shock  front,  at  77  =  t)s ,  are  the  same  as  in  Ref.  13: 


m,) 


2ns -r+i 


y  +  1 


(A8) 
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S(ils)  = 


2/07,  +i) 

7  +  1 


(A9) 


The  self-similar  coordinate  of  the  shock  front,  77  s ,  is  not  known  in  advance.  The  two 
constants,  ps  and  S0  in  Eq.  (A7),  must  be  chosen  to  select  a  solution  of  Eqs.  (A3),  (A4)  that 


satisfies  two  boundary  conditions  (A8),  (A9).  This  is  done  by  the  shooting  method.  For  instance, 
one  can  select  a  trial  value  of  S0  and  use  the  asymptotic  expressions  (A6)-(A7)  to  start  integration 
of  Eqs.  (A3),  (A4)  from  some  small  but  finite  77  until  U ( 77 )  equals  (2t]  -  7  + 1 )/( 7  + 1 )  so  that  the 
boundary  condition  (A8)  is  satisfied  at  some  value  of  77  approximating  J)s .  Then  we  see  if  the  left- 
hand  side  of  Eq.  (A9)  is  larger  or  smaller  than  its  right-hand  side,  adjust  the  value  of  Sf) 

accordingly,  and  the  procedure  is  repeated  until  it  converges.  This  method  was  used  to  plot  the  self¬ 
similar  profiles  shown  in  Figs.  1,  2. 

Appendix  II.  Boundary  conditions  at  the  ablation  front. 

We  derive  the  boundary  conditions  at  the  ablation  front  for  the  case  of  time-dependent 
acceleration  and  ablation  velocity,  from  those  given  in  Ref.  9.  Our  notation  is  slightly  different  from 
that  of  Ref.  9.  We  suppose  the  plasma  to  be  driven  in  the  positive  x  direction,  and  the  small 
perturbations  °c  exp  (iky)  (the  exponential  factor  is  omitted  below)  to  be  imposed  in  the  x-y  plane. 

In  the  reference  frame  of  the  ablation  front,  velocity  of  the  plasma  is  thus  negative,  equal  to 
-  v,  and  ~v 2  upstream  and  downstream  from  the  ablation  front,  respectively.  By  definition, 
v,  =  va ;  due  to  mass  conservation  in  the  quasi-steady  ablation  front,  p2  /  pl  =  v1  /  v2  =  rD  « 1 .  The 
physical  assumptions  of  Ref.  9  remain  unchanged.  Both  characteristic  velocities,  v,  andv2,  are 
much  less  than  the  respective  sound  speeds;  thermal  conduction  is  strong  downstream  and 
negligible  upstream,  the  ablation  front  is  an  isotherm.  Then  the  linearized  equations  of  continuity, 
motion,  and  energy  yield  the  following  boundary  conditions: 


(Pi  =  PA',  1  -  Pi^xi  > 

at 

8vyi  -  Svy2  =  -ikv2(l  -  rD  )8xa , 

<S/h  -  Sp2  -  2 PjVj  (8vxl  -  8vxl )  -  (p,  -  p2)g8xa  =  0, 


(A10) 

(All) 
(A  12) 


8vxl-8vx2  =  -kv2(l-rD)8xa 


(A13) 


Here,  8p,8vx  and<5vv  are  perturbations  of  pressure,  axial  and  lateral  velocity,  and  the  subscripts  1 

and  2  refer  to  the  high-density  (upstream)  and  low-density  (downstream)  plasmas  separated  by  the 
ablation  front,  respectively.  Under  the  assumptions  outlined  above,  density  perturbations  drop  out 
from  the  boundary  conditions. 

Substituting  (A13)  into  (A10),  we  obtain  Eq.  (9).  From  (A12),  (A13),  we  find: 


Pi 


Pi 


(A  14) 


We  have  to  exclude  the  variables  referring  to  the  downstream  state.  Assuming  that  no 
perturbations  are  incident  onto  the  ablation  front  from  the  plasma  corona  (from  downstream),  one 
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can  find  the  dependence  between  the  downstream  perturbation  amplitudes 


Svx2,  and  8v  y2 


from  the  dispersion  relations  for  the  outgoing  waves.  After  some  algebra,  we 


obtain: 


kSp 


d 


2  =  —rn  — 8vr ,  -ikv  8v 


a'~'v  y2  ■ 


(A  15) 


Pi  dt 

In  (A15),  we  express  dvx2  from  (A13),  8vy2  from  (All),  substitute  the  resulting  expression  for 
8}2  into  (A14),  and  simplify  it  with  the  aid  of  (9).  We  are  left  with  a  boundary  condition  that  relates 
only  the  perturbations  of  the  upstream  flow  variables,  so  that  the  subscript  1  in  the  perturbation 
amplitudes  could  be  omitted: 


k8p 

~pT 


+  (1  ~rD) 


(1  0>)  I  2  2  7 

- k  va  -  gk 


5xa  +  kva(8vx+i8vv)  +  rl 


PC  d  ^ 

»  5*.-* 

\  J 


-0.  (A  16) 


This  is  a  more  general  form  of  the  boundary  condition  sought  for.  It  could  be  used  to  describe  the 
effect  of  perturbations  incident  to  the  ablation  front  from  the  dense  plasma  (in  this  case, 

8p,8vx,  and  Sv  are  independent  variables).  The  stabilizing  “rocket  effect”  of  Ref.  9  is  expressed  by 

the  term  proportional  to  k2v2  /  rD  in  the  left-hand  side  of  Eq.  (A16):  a  displacement  of  the 
ablation  front  generates  a  negative  contribution  to  the  pressure  response  x  ~dxa  that  tends  to 


decrease  the  displacement. 

To  derive  Eq.  (10),  we  have  to  make  an  additional  assumption  that  no  perturbations  are 

incident  onto  the  ablation  front  from  the  upstream,  thus  neglecting  the  vorticity  perturbations  that 
arrive  with  the  plasma  particles.  Then  the  perturbation  amplitudes  8p,  8vx,  and  8v v  are  connected 


via  the  dispersion  relation  for  sonic  waves  propagating  upstream: 

—=irSv*+kvaSv*’  Sv--i8vx.  (A17) 

Pi  dt 

Substituting  (A17)  into  (A16),  we  obtain  (10). 


16 


References 


I  S.  E.  Bodner,  D.  G.  Colombant,  J.  H.  Gardner,  R.  H.  Lemberg,  S.  P.  Obenschain, 

L.  Phillips,  A.  J.  Schmitt,  J.  D.  Sethian,  R.  L.  McCrory,  W.  Seka,  C.  P.  Verdon, 

J.  P.  Knauer,  B.  B.  Afeyan,  and  H.  Powell,  Phys.  Plasmas  5,  1901  (1998). 

K.  O.  Mikaelian,  Phys.  Rev.  Lett.  48,  1793  (1982) 

N.  A.  Inogamov,  Prikl.  Mekh.  Tekh.  Liz.  (in  Russian)  5,  110  (1985). 

4  N.  A.  Inogamov,  JETP  84,  746  (1997);  Phys.  Rev.  E  57,  3352  (1998) 

5  R.  Ishizaki  and  K.  Nishihara,  Phys.  Rev.  Lett.  78,  1920  (1997). 

6  R.  J.  Taylor,  A.  L.  Velikovich,  J.  P.  Dahlburg,  and  J.  H.  Gardner,  Phys.  Rev.  Lett.  79, 
1861  (1997);  A.  L.  Velikovich,  J.  P.  Dahlburg,  J.  H.  Gardner,  and  R.  J.  Taylor,  Phys. 
Plasmas  5,  1491  (1998). 

7V.N.  Goncharov,  R.  Betti,  R.L.  McCrory,  C.  Cherfils,  Bull.  Am.  Phys.  Soc.  43,  1667 
(1998);  V.  N.  Goncharov,  Phys.  Rev.  Lett.,  to  be  published  (1999). 

8  J.  Sanz,  Phys.  Rev.  Lett.  73,  2700  (1994). 

9  A.  R.  Piriz,  J.  Sanz  and  L.  L.  Ibanez,  Phys.  Plasmas  4,  1117  (1997). 

10  A.  L.  Velikovich,  L.  L.  Cochran,  and  J.  Davis,  Phys.  Rev.  Lett.  77,  853  (1996); 

A.  L.Velikovich,  L.  L.  Cochran,  J.  Davis,  and  Y.  Chong,  Phys.  Plasmas  5,  3377 
(1998). 

II  J.  H.  Hammer,  J.  L.  Eddleman,  P.  Springer,  M.  Tabak,  A.  Toor,  K.  Wong, 

G.  B.  Zimmerman,  C.  Deeney,  R.  Humphreys,  T.  J.  Nash,  T.  W.  L.  Sanford, 

R.  B.  Spielman  and  J.  S.  DeGroot,  Phys.  Plasmas  3,  2063  (1996). 

12  N.  A.  Krasheninnikova,  Izv.  Akad.  Nauk  SSSR,  Otd.  Tekh.  Nauk  8,  22  (1955)  (in 
Russian);  a  brief  description  of  this  solution  and  references  to  subsequent  Russian 
publications  are  given  by  L.  I.  Sedov,  Similarity  and  Dimensional  Methods  in 
Mechanics  (CRC  Press,  Boca  Raton,  1993),  Chap.  IV,  Sec.  6. 

13  S.  M.  Gol’berg  and  A.  L.  Velikovich,  Phys.  Lluids  B  5,  1164  (1993). 

14  R.  L.  McCrory,  L.  Montierth,  R.  L.  Morse,  and  C.  P.  Verdon,  in  Laser  Interaction  and 
Related  Plasma  Phenomena  (Plenum,  New  York,  1981),  Vol.  5,  pp.  713-742;  R.  Betti, 
R.  L.  McCrory,  and  C.  P.  Verdon,  Phys.  Rev.  Lett.  71,  3131  (1993). 

15  L.  J.  Slater,  “Confluent  Hypergeometric  Lunctions”,  in  Handbook  of  Mathematic cd 
Functions ,  ed.  by  M.  Abramowitz  and  I.  A.  Stegun  (Dover,  New  York,  1972),  pp.  503- 
536. 

16  V.  N.  Goncharov,  R.  Betti,  R.  L.  McCrory,  P.  Sorotokin,  and  C.  P.  Verdon,  Phys. 
Plasmas  3,  1402  (1996). 

17  S.  E.  Bodner,  Phys.  Rev.  Lett.  33,  761  (1974). 

18  N.  C.  Lreeman,  Proc.  Roy.  Soc.  A  228,  341  (1955);  P.  M.  Zaidel,  J.  Appl.  Math.  Mech. 
24,  316  (1960). 


17 


19  M.  H.  Emery,  J.  H.  Gardner,  and  J.  P.  Boris,  Appl.  Phys.  Lett.  41,  808  (1982). 

20 

E.  S.  Oran  and  J.  P.  Boris,  Numerical  Simulations  of  Reactive  Flow  (Elsevier,  New 
York,  1987),  and  references  therein. 

21  F.  F.  Grinstein,  and  C.  R.  DeVore,  Phys.  Fluids  8,  1237  (1996). 

2“  A.  M.  Landsberg,  J.  P.  Boris,  W.  Sandberg,  and  T.  R.  Young,  Jr.,  AIAA  Paper  95- 
0047,  American  Institute  of  Aeronautics  and  Astronautics,  Reston,  VA  20101-4344, 
January  1995. 

23  J.  H.  Gardner,  S.  E.  Bodner,  and  J.  P.  Dahlburg,  Phys.  Fluids  B  3,  1070  (1991). 

24  R.  J.  Taylor,  J.  P.  Dahlburg,  A.  Iwase,  J.  H.  Gardner,  D.  E.  Fyfe,  and 

O.  Willi,  Phys.  Rev.  Lett.  76,  1643  (1996);  C.  J.  Pawley,  K.  Gerber,  R.  H.  Lehmberg, 
E.  A.  McLean,  A.  N.  Mostovych,  S.  P.  Obenschain,  J.  D.  Sethian,  V.  Serlin, 

J.  A.  Stamper,  C.  A.  Sullivan,  S.  E.  Bodner,  D.  Colombant,  J.  P.  Dahlburg, 

A.  J.  Schmitt,  J.  H.  Gardner,  C.  Brown,  J.  F.  Seely,  T.  Lehecka,  Y.  Aglitskiy, 

A.  V.  Deniz,  Y.  Chan,  N.  Metzler,  and  M.  Klapisch,  Phys.  Plasmas  4,  1969  (1997). 

25  J.  P.  Boris  and  D.  L.  Book,  J.  Comput.  Phys.  11,  38  (1973),  also  see  “Solution  of  the 
Continuity  Equation  by  the  Method  of  Flux-Corrected  Transport,”  Methods  in 
Computational  Physics ,  v.  16  (Academic  Press,  New  York,  1976),  pp.  85-129. 

26  L.  Spitzer  and  R.  Harm,  Phys.  Rev.  89,  977  (1953). 

27 

S.  E.  Bodner,  D.  Colombant,  K.  Gerber,  R.  H.  Lehmberg,  E.  A.  McLean, 

A.  N.  Mostovych,  M.  S.  Pronko,  C.  J.  Pawley,  A.  J.  Schmitt,  J.  D.  Sethian,  V.  Serlin, 
J.  A.  Stamper,  C.  A.  Sullivan,  J.  P.  Dahlburg,  J.  H.  Gardner,  Y.  Chan,  A.  V.  Deniz, 

J.  Hardgrove,  T.  Lehecka,  and  M.  Klapisch,  Phys.  Plasmas  3,  2098  (1996). 


18 


0.0  0.5  1 .0 


Fig.  1.  Self-similar  profiles  of  density  N  (thick  lines),  pressure  P  and  temperature  0  for  the 
shock-piston  flows  driven  by  a  linearly  rising  pressure  (the  self-similar  coordinate  77  is  normalized 
with  respect  to  the  shock  coordinate,  rjs ):  (a)  tailored  (parabolic)  initial  density  profile  (x  =  2), 

deceleration  (m  =  3/4);  (b)  uniform  density  profile  (%  =  0),  acceleration  (m  =  3/2). 
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0.0  0.5  1.0 


Fig.  2.  Self-similar  density  profile  of  Fig.  1(a)  (thick  solid  curve)  compared  to  the  density  profiles 
obtained  in  a  simulated  ablatively  driven  shock-piston  flow  that  approximately  corresponds  to  the 
same  parameters  =  2,  m  =  3/4).  Simulation  results  (13  thin  solid  and  dotted  curves)  are 
calculated  during  the  shock  transit  time  for  this  target,  between  1.5  and  3.5  ns.  The  flow  profiles 
are  seen  to  be  approximately  self-similar. 
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kL(t) 


kL(t) 


Fig.  3.  Evolution  of  perturbations  of  the  ablation  front  given  by  (13)  for  rD  =  0.05,  A  =  0.9 ,  m 
3/4  and  m  =  3/2  (thin  and  thick  solid  curves,  respectively):  (a)  no  mass  ablation,  v  ♦  0;  (b)  with 
mass  ablation,  v=  0.05. 


21 


l[TW/cm?]  p[g/crrf] 
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Fig.  4.  (a)  Initial  density  profiles  in  targets:  parabolic  (1)  and  low-density  uniform  (2)  profiles 
followed  by  a  solid  payload;  all-uniform  low-density  (3)  and  solid  plastic  density  (4)  profiles,  (b) 
Laser  pulse  shape  selected  to  make  ablative  pressure  increase  linearly  with  time  from  f,  =  0.2  ns  to 

t2  =  4  ns  :  p  /“ ,  where  a  =  0.75,  so  we  choose  1  (t)  =  - t)  +  /“  (t  —  tl)~\(t2  -  tx)\a , 

where  /,  =  3  ■  101 2  W/cm 2  and  =  1 .9  ■  101 3  W/cm 2 .  The  pressure  rise  is  indeed  approximately 

linear  in  time. 
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P[Mbar] 


Fig.  5.  x-t  diagrams  for  plasma  particles  labeled  by  their  initial  mass  coordinates  for:  (a)  target  1 
with  a  parabolic  density  profile;  (b)  target  2  with  a  low-density  uniform  profile;  both  followed  by  a 
solid  density  payload.  The  upper  envelope  curve  of  the  high-density  area  is  the  ablation  front.  At 
early  time,  in  the  parabolic  density  target,  this  curve  is  convex,  which  implies  deceleration. 
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Increasing  Laser  power 
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Fig.  6.  Time  evolution  of  areal  mass  variation  ((im^g/cm3)  in  targets  1  to  4  driven  by  the  laser 
pulse  of  Fig.  4(b),  with  a  0.2%  rms  single-mode  lateral  intensity  variation  at  X  -  15  ftm. 
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Increasing  Laser  power 
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Fig.  7.  Time  evolution  of  areal  mass  variation  ((im^g/cm3)  in  targets  1  to  4  driven  by  the  laser 
pulse  of  Fig.  4(b),  with  a  0.2%  rms  single-mode  lateral  intensity  variation  at  X  =  30  (im. 
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Increasing  Laser  power 
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Fig.  8.  Time  evolution  of  areal  mass  variation  (jim^g/cm3)  in  targets  1  to  4  driven  by  the  laser 
pulse  of  Fig.  4(b),  with  a  0.2%  rms  single-mode  lateral  intensity  variation  at  A  =  60  Fm- 
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Fig.  9.  Time  evolution  of  areal  mass  variation  ((im^g/cm3)  in  targets  1  to  4  driven  by  the  laser 
pulse  of  constant  intensity,  with  a  0.2%  rms  single-mode  lateral  intensity  variation  at  A  =  15  ftm . 
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Fig.  10.  Time  evolution  of  areal  mass  variation  (pm^g/cm3)  in  targets  1  to  4  driven  by  the  laser 
pulse  of  constant  intensity,  with  a  0.2%  rms  single-mode  lateral  intensity  variation  at  A  =  30  Tm- 
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Constant  Laser  Power 
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Fig.  11.  Time  evolution  of  areal  mass  variation  (pm^g/cm3)  in  targets  1  to  4  driven  by  the  laser 
pulse  of  constant  intensity,  with  a  0.2%  rms  single-mode  lateral  intensity  variation  at  X  =  6C  pm. 
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Parabolic  Density;  Laser  Perturbation 
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Fig.  12.  Density  perturbation  maps  of  p  -  H  in  a  tailored  density  target  1  driven  by  a  laser 


pulse  of  Fig.  4(b)  with  a  15  (im  perturbation  wavelength  at  different  times:  (a)  1.5  ns;  (b)  2  ns;  (c) 
4.5  ns;  (d)  10  ns. 


